Astron. Nachr. / AN , No. , 789-17921(2010) / DOI please setDOI! 



Eclipsing binaries with pulsating components: CoRoT 102918586 
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We present the preliminary results of the study of an interesting target in the first CoRoT exo-planet field (IRal): 
CoRoT 102918586. Its light curve presents additional variability on the top of the eclipses, whose pattern suggests multi- 
frequency pulsations. The high accuracy CoRoT light curve was analyzed by applying an iterative scheme, devised to 
disentangle the effect of eclipses from the oscillatory pattern. In addition to the CoRoT photometry we obtained low reso- 
lution spectroscopy with the AAOmega multi-fiber facility at the Anglo Australian Observatory, which yielded a spectral 
classification as F0 V and allowed us to infer a value of the primary star effective temperature. The Fourier analysis of the 
residuals, after subtraction of the binary light curve, gave 35 clear frequencies. The highest amplitude frequency, of 1.22 
c/d, is in the expected range for both 7 Dor and SPB pulsators, but the spectral classification favors the first hypothesis. 
Apart from a few multiples of the orbital period, most frequencies can be interpreted as rotational splitting of the main 
frequency (an 1 = 2 mode) and of its overtones. 
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1 Pulsators in CoRoT eclipsing binaries 

The space mission CoRoT0has already provided high qual- 
ity photometry and continuous monitoring of about 140000 
stars. Its photometric accuracy ranges from a few to a few 
hundred parts per million (Auvergne et al. 2009), depend- 
ing on star magnitude and on the observing channel (Aster- 
oseismo or Exoplanet); the duration of the observing runs 
varies from 20 to 150 days (for short or long runs), with a 
typical duty cycle > 90%. With these figures there is no sur- 
prise that most CoRoT targets show some kind of variabil- 
ity, often of quasi-periodic character. In particular CoRoT 
has discovered and monitored hundreds of eclipsing bina- 
ries (EBs). Thanks to the Corot Variability Classifier (CVC, 
Debosscher et al. 2009), and a further detailed screening, 
we selected a group of EBs whose light curves present other 
periodic phenomena besides eclipses, possibly due to a pul- 
sating component. Even though the analysis is made more 
complex by the superimposed phenomena, the benefit is the 
access to independent, and somewhat complementary, sour- 
ces of information on the stellar fundamental parameters 
(double-lined spectroscopic binaries provide, for instance, a 
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pure geometrical determination of masses and radii, which 
are most useful in asteroseismic studies). 

CoRoT 102918586 is one of the most promising targets 
of our sample, because of its brightness (V~ 12.4), which 
allows spectroscopic follow-up, and of its light curve shape 
clearly showing modulated oscillations and narrow eclipses 

(Fig.LD. 

2 Method of analysis 

All light curves of our sample need a preliminary treat- 
ment to remove possible instrumental trends and obvious 
outliers. When none of the two periodic phenomena (pul- 
sation and orbit) can be treated as a perturbation, as in the 
case of CoRoT 102918586, we apply the following iterative 
scheme: 

A. Fourier analysis (with Period04, Lenz & Breger 2005), 
of a time series obtained after removal of the eclipses. 
The frequency significance empirical criterion of Breger 
etal. (1993) is used. 

B. Prewhitening of the de-trended full light curve (here- 
after LC0) with the derived pulsation frequencies, pro- 
viding a first-approximation binary light curve (LCB1) 

C. Light curve solution based on LCB 1 (or a conveniently 
resampled version) with WD-PHOEBE (Prsa & Zwit- 
ter 2005) and computation of the residuals from the full 
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Fig. 1 The de-trended re-binned to 512 s (for better readability) light curve of CoRoT 102918586 and the final fit (full 
line) obtained by summing the binary model and the harmonic signal from the 35 detected frequencies. 



light curve LC0), yielding a new (and continuous) time 
series to submit to Period04. 
D. Second stage of prewhitening of LC0, providing LCB2, 
on which the PHOEBE solution can be refined. Steps 
B and C need to be repeated at least once, to refine the 
PHOEBE solution, but the loop can be repeated as many 
times as necessary. 

3 Observations of CoRoT 102918586 
3.1 CoRoT photometry 

CoRoT 102918586 was observed for about 57 days (from 
HJD = 2454135.056 to 2454192.499) during CoRoT's ini- 
tial run in the exoplanet field (IRal), with a time sampling 
of 512 s for the first 27 days and then of 32 s for the rest of 
the record. Little information was known about this target 
prior to the CoRoT mission, and its binary nature was re- 
vealed by the CoRoT Variability Classifier (Debosscher et 
al. 2009). 

The light curve has a total of 87355 points, and is re- 
corded in "chromatic" light, i.e. in three uncalibrated bands 
unevenly spanning the wavelength interval 400-900 nm, ob- 
tained by the insertion of a bi-prism just before the CCD 
of the exo-planet channel. In this particular case, the three 
chromatic curves do not present any notable differences and 
we added them in order to obtain a single light curve with 
increased Signal to Noise Ratio (SNR). After the cleaning of 



obvious outliers, the removal of a slow decreasing trend due 
to some instrumental decay and homogenization of the sam- 
pling, we obtained the normalized curve shown in Fig. [T] 
The examination of the contaminants in the mask used to in- 
tegrate the flux reveals that the brightest one is 3 magnitudes 
fainter than the target, and is, therefore, negligible. The light 
curve displays periodic variations in intensity presenting a 
beating pattern, with a maximum peak to peak amplitude 
of about 25 mmag that could be pulsations of one or both 
the components. Moreover, the eclipsing binary nature can 
clearly be seen thanks to deep regularly spaced minima. 

3.2 Ground-based additional information 

The ExoDat databas^l provides for this target the color in- 
formation collected in table Q] In addition, we took advan- 

Table 1 Colour photometry for CoRoT 102918586 





Harris B 


Harris V 


Sloan r' 


Mag 


12.842 ± 0.03 


12.425 ±0.029 


12.244 ± 0.048 



tage of the AAOmega multi-object facilities at the Anglo- 
Australian Observatory to acquire a first spectroscopic snap- 



The Exo-Dat database is operated at LAM-OAMP, Marseille, France, 
on behalf of the CoRoT/Exoplanet program 
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shot of the CoRoT fields. Low resolution spectroscopy (with 
R « 1300) of the CoRoT targets in the Anticenter fields 
was carried out with the multi-fiber spectrograph during two 
observing runs, in January 2008 and December 2008 - Jan- 
uary 2009. The spectrum obtained for CoRoT 102918586 is 
shown in (Fig. [2j together with its best fit F0 V template. 

4 The eclipsing binary model 

Using the Phase Dispersion Minimization method (Stelling- 
werf, 1978) and refining the result with PHOEBE, we de- 
rived the orbital period P = 8'?78248. Although a period 
of a half cannot be excluded from photometry alone, con- 
sidering that: i) there is no clear evidence of a secondary 
minimum corresponding to a 4*? 39 period ; ii) after itera- 
tion (see section |2), a slight difference can be seen between 
minima (in the sense of larger dispersion in the secondary 
eclipse) in accordance with the 81*78 period, we decided to 
keep the value P = 8^78. A final answer will come from 
the high resolution spectroscopy we plan to acquire for this 
target. 

After pre-whitening with the fifteen frequencies detec- 
ted in step A of the iterative procedure (see section [2j, we 
proceeded to the first LC solution. The original curve was 
re-sampled to 512 s bins in the out-of-eclipse intervals (to 
decrease computing time and the numerical weight of the 
large number of points in the flat sections). Step D (2nd 
pre-whitening) was then made. Degeneracy cannot be avoid 
while solving the inverse problem for eclipsing binaries, all 
the more so when there is only one input light curve, and 
the value of some parameters need to be assumed. Here we 
set the primary temperature T e s,i at the one deduced from 
spectroscopy (± 200 K) and assumed ms/mi = 1 (it can- 
not be derived from the LC solution of detached systems). 
Albedo and gravity darkening were fixed at the theoreti- 
cal values for radiative envelopes (1.0), and accordingly we 
used a logarithmic limb darkening law. The choice of the pe- 
riod implied a circular orbit (equidistant minima). The final 
PHOEBE solution is shown in Fig. [3] and the main param- 
eters are given in Table [2] (where r is the fractional radius). 
The global fit (binary + pulsations) is shown in Fig. Q] All 
fits are obtained with the updated version of PHOEBE (Jan 
2010) including model atmospheres and limb darkening ta- 
bles in the CoRoT passband (Maceroni et al. 2009). 

Table 2 Main parameters of the binary model. All errors 
are formal errors derived from the least-squares fit. 



System 



Primary 



Secondary 



Teff (K) 

L/ Ltot 
r 



86.053 ± 0.001 



7300" 
0.5283 ±0.0001 
0.04081 ± 0.00001 



7350 ± 10 
0.4717 ±0.0001 
0.0381 ±0.0001 



a) Fixed value 
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Fig. 2 Low resolution spectrum (R^1300) of CoRoT 
102918586 (black line), observed with the AAOmega multi- 
object facility at the Anglo-Australian Observatory. The 
best fitting F0 V template has been superimposed with a 
broken line (T off = 7300 K, log g=4.0). 
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Fig. 3 The residuals pre-whitened with the final 35 fre- 
quencies with S/N > 4, and yielding the eclipsing binary 
light curve (LCB2). The original ^75000 points are reduced 
in the LC solution to ^21000, by resampling to 512 s bins 
outside eclipses. The final solution is superimposed with a 
continuous line. 



5 Pulsation analysis 

The Fourier analysis of the residuals, after subtraction of 
LCB2 fit (Fig. [3]), yields 35 significant frequencies. The first 
thirteen are shown in Table [3] and Fig. [4] The frequency 
range, around 1 c/d, and the spectral classification suggest 
the presence of a 7 Dor component. The eclipses being par- 
tial, we cannot tell which one is the pulsating component 
(though a decrease of dispersion at phase zero favors the 
hypothesis of a pulsating primary). 

The straightforward interpretation of the frequency pat- 
tern is that we are observing the rotational splitting of an 
1=2 g-mode around F3 and its overtones (2i<3 and 3F3), and 
a few multiples of the orbital frequency F or b = 0.1139 c/d, 
which could be related to small systematic deviations of the 
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Table 3 The first thirteen frequencies, ordered by am- 
plitude (as detected by Period04 in the normalized light 
curve prewithened with the binary signal), their correspond- 
ing phases and S/N ratios. The numbers in brackets are the 
errors on the last digit. The S/N ratio was computed over an 
interval of 5 c/d. 





F (c/d) 


Ampl -10 a 


Phase (2tv rad) 


S/N 


Fx 


1.22452 (2) 


4.30116 (1) 


0.1576 (3) 


24.7 


F-2 


1.12666 (3) 


3.05882 (1) 


0.1223 (4) 


21.7 


F 3 


1.17256 (3) 


2.64945 (1) 


0.6155 (5) 


24.2 


F 4 


0.94653 (4) 


1.43687 (1) 


0.2809 (9) 


16.5 


ft 


2.3512(1) 


0.72474 (1) 


0.101 (2) 


11.8 


F 6 


0.4572 (1) 


0.66486 (1) 


0.432 (2) 


8.51 


F 7 


0.2277 (1) 


0.53671 (1) 


0.286 (2) 


7.00 


ft 


2.39622 (1) 


0.49724 (1) 


0.903 (3) 


9.77 


ft 


0.0519 (2) 


0.50943 (1) 


0.251 (3) 


6.72 


Fxo 


2.4490 (2) 


0.41644 (1) 


0.097 (3) 


9.65 


Fxx 


0.6833 (2) 


0.44263 (1) 


0.055 (3) 


7.96 


F 12 


2.3009 (2) 


0.42179 (1) 


0.944 (3) 


9.33 


F 13 


0.9118 (2) 


0.32531 (1) 


0.901 (4) 


6.49 



Table 4 Frequency interpretation. 





F (c/d) 


Comment 




F (c/d) 


Comment 


F25 


1.074 


ql, m=-2 


F24 


2.2499 


q2, m=-2 


ft 


1.12666 


ql, m=-l 


FX2 


2.3009 


q2, m=-l 


ft 


1.17256 


ql, m=0 


F 5 


2.3512 


q2, m=0 


ft 


1.22452 


ql, m=+l 


ft 


2.39622 


q2, m=+l 


fto 


1.2756 


ql, m=+2 


fto 


2.4490 


q2, m=+2 


fto 


3.476 


q3, m=-l 


F 7 


0.2277 


2 -Porb 


F23 


3.5254 


q3, m=0 


ft 


0.4572 


4F orb 


F27 


3.5740 


q3, m=+l 


Fxx 


0.6833 


6 ftrb 


F 35 


3.622 


q3, m=+2 


F 13 


0.9118 


8 -Forb 



binary light curve fit (Table|4]i. The mean frequency splitting 
of about 0.05 c/d corresponds approximatively to 0. Sorb- 
in a synchronized system - as a circular orbit suggests - F or b 
would also be the frequency of rotation. The splitting is in- 
deed of the right order to be of rotational origin. Note that 
Fz/Foth ?» 10, a similar ratio for tidally excited modes be- 
ing already found in the SPB HD 177863 (Willems & Aerts 
2002) and HD 174884 (Maceroni et al. 2009) 

We expect to see I = 2 modes excited in presence of 
a dynamical tide or an equilibrium tide, or both (e.g. Zahn, 
2005). However, the first mechanism requires an eccentric 
orbit, whereas the second requires tidally deformed compo- 
nents. None of these conditions is fulfilled by the present 
binary model. An alternative model could be built in the hy- 
pothesis of an orbital period half the assumed value and a 
particular orientation in space of an eccentric orbit, so that 
only one eclipse would occur in a system with similar com- 
ponents (SB2). A similar configuration is at the origin, for 
instance, of the light curve of HD 174884 (Maceroni et al. 
2009) where twin components produce minima with a ratio 
of depths of about a hundred. 




12 3 4 

Frequency (c/d) 



Fig. 4 Fourier spectrum before pre-whitening, with indi- 
cations of the first thirteen frequencies of Table [5J and the 
spectral window. 

6 Conclusions 

The preliminary results on C0R0T 102918586 are based on 
a single pass-band (though excellent) photometry and low 
resolution spectra, nevertheless the combination of pulsa- 
tion and binarity yields already a rich harvest of informa- 
tion. Obviously we need high resolution and time resolved 
spectroscopy to definitely solve the remaining ambiguities 
of this interesting system. This target is included with high 
priority in our spectroscopic follow-up programs. 

Acknowledgements. The research leading to these results has re- 
ceived funding from the Italian Space Agency (ASI) under con- 
tract ASI/INAF 1/015/07/00 in the frame of the ASI-ESS project, 
from the European Research Council under the European Com- 
munity's Seventh Framework Programme (FP7/2007-2013)/ERC 
grant agreement n° 227224 (PROSPERITY), as well as from the 
Research Council of K.U.Leuven grant agreement GOA/2008/04. 
CM acknowledges the European Helio and Asteroseismology Net- 
work (HELAS), a major international collaboration funded by the 
European Commission's Sixth Framework Programme, for sup- 
porting her contribution to the conference. 

References 

Auvergne, M., Bodin, P., Boisnard, L. et al.: 2009, A&A 506, 41 1 
Breger, M., Stich, J., Garrido, R. et al.: 1993, A&A, 271, 482 
Damiani, C, Maceroni, C, Cardini, D. et al.: 2010, Ap&SS (Helas 

3 Conf. Proc.) in press 
Debosscher J., Sarro, L. M., Lopez, M. et al.: 2009, A&A, 506, 

519 

Lenz P., Breger M.: 2005, CoAst, 146, 53 

Maceroni C, Montalban, J., Michel, E. et al.: 2009 A&A, 508, 
1375 

Prsa, A., Zwitter, T: 2005, ApJ, 628, 426 

Stellingwerf, R. F. : 1978, ApJ, 224, 953 

Willems, B., Aerts, C: 2002, A&A, 384, 441 

Zahn J-P: 2005, in "Tidal Evolution and Oscillations in Binary 

Stars: Third Granada Workshop on Stellar Structure", ASPC 

333,4 



© 2010 WILEY- VCH Verlag GmbH & Co. KGaA, Wcinhcim 



www.an-journal.org 



